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ABSTRACT 

We construct new galaxy angular power spectra Ci based on the extended, updated 
and final SDSS II Luminous Red Galaxy (LRG) photometric redshift survey-MegaZ 
(DR7). Encapsulating 7746 deg 2 we utilise 723,556 photometrically determined LRGs 
between 0.45 < z < 0.65 in a 3.3 (Gpc h^ 1 ) 3 spherical harmonic analysis of the 
galaxy distribution. By combining four photometric redshift bins we find preliminary 
parameter constraints of fb = flb/D, m = 0.173 ± 0.046 and fl m = 0.260 ± 0.035 as- 
suming Hq = 75 km s _1 Mpc -1 , n s = 1 and ftk = 0. These limits are consistent 
with the CMB and the previous data release (DR4). The Cg are sensitive to redshift 
space distortions and therefore we also recast our constraints into a measurement of 
(3 ss S!" n 55 /6 in different redshift shells. The robustness of these power spectra with 
respect to a number of potential systematics such as extinction, photometric redshift 
and ANNz training set extrapolation are examined. The latter includes a cosmologi- 
cal comparison of available photometric redshift estimation codes where we find ex- 
cellent agreement between template and empirical estimation methods. MegaZ DR7 
represents a methodological prototype to next generation surveys such as the Dark 
Energy Survey (DES) and, furthermore, is a photometric precursor to the spectro- 
scopic BOSS survey. Our galaxy catalogue and all power spectra data can be found 
at |http://zuserver2. star. ucl.ac.uk/~sat/ MegaZ /McgaZD R7.tar.gz| 

Key words: Large Scale Structure, Galaxy Clustering, Dark Energy, Cosmological 
Parameters and Photometric Redshifts. 



1 INTRODUCTION 

The analysis of the statistical distribution of fluctuations in 
the Universe is a potent method for constraining theories or 
components within Cosmology. In fact, the power spectrum 
will fully describe these variations, which are predicted by 
theory, if they are given by a Gaussian random field. The 
Cosmic Microwave Background (CMB) has been a great ex- 
ample of this principle in action with recent high precision 
measurements ( Komatsu et al.|2010 1 confirming that a clear 
and consistent picture of cosmology is emerging. It is desir- 
able however to test this picture with additional and inde- 
pendent data that explores a contrasting epoch of cosmic 
evolution and breaks the parameter degeneracies that exist 
from a single probe of the early Universe. A galaxy redshift 



survey is therefore a powerful tool in Cosmology (Peebles 



19731. In addition, this late-time galaxy distribution is sen- 



sitive to the emergence of dark energy ( |Riess et al.||1998| 



of structure enables a test of gravity ( 


Jain & Zhang 2007 


|Huterer & Linder||2007[ |Thomas et al. 


2009') and the mass 


of the neutrino (Hu et al.||1998 and Thomas et al.|2010l. 



The structure and aim of this paper is as follows: To 
construct and present the angular power spectra Ce of the 
new SDSS Luminous Red Galaxy (LRG) photometric sur- 
vey, along with the associated error and individual cosmolog- 
ical constraints. Specifically, we determine the colour, red- 
shift and angular selection functions that define the survey 
in Section [2] The spherical harmonic analysis is described 
in Section [3] and Section [4] The cosmological constraints 
inferred and the potential systematics of the data set are 
discussed in Section [5] and Section [6j respectively. This last 
section also includes a cosmological comparison of different 
photometric redshift methods. Extended and combined cos- 
mological implications are to be presented in a companion 
paper with the likelihood. 



2 DATA 

The development of galaxy surveys over the past few years 
reflects the balance between observational technology and 
gains in cosmological parameter estimation. This has at 
present culminated in the impressive 2-degree Field Galaxy 



2 Thomas et al. 



Redshift Survey (2dFGRS - [Colless et aT^|200T} and the 



Sloan Digital Sky Survey (SPSS - |Yorket al.||2000[ ). How- 
ever, the acquisition of a vast number of precise redshifts 
through spectroscopy is an expensive, challenging and time 
consuming task. An alternative method is to use photomet- 



ric redshifts (E.g. Csabai et al. 2003 I resulting from observa- 
tions of broadband galaxy colours through a series of filters. 
The motivation is that a decrease in redshift accuracy is out- 
weighed by measurements of a vast number of galaxies over 
a wide area of the sky, therefore encompassing a large cos- 
mic volume. Photometric redshift surveys have been shown 
to be competitive by Blake et al. 2007] and |Padmanabhan| 
|et al.|[2007| and here we follow these papers. Upcoming sur- 



veys, such as the Dark Energy Survey (The Dark Energy 
Survey Collaboration 20051, are heavily based on this effi- 



ciency principle. 

We therefore aim to analyse the clustering of the lat- 
est and final SDSS II photometry given by Data Release 7 
(DR7) - |Abazajian et al.| j p009t ) . The « 1.5 million LRG cat- 
alogue (MegaZ-LRG DR7) is produced as an updated ver- 
sion of MegaZ-LRG ( |Collister et al.|2007[ ). These LRGs are 
old red elliptical galaxies that provide a clean and consistent 
galaxy sample. With a stable spectral energy distribution 
(SED) and a sharp 4000A break they provide good photo- 
metric redshift estimates. Furthermore, they are known to 
strongly trace the underlying mass density; a distribution 
we are striving to quantify. Also, being among the brightest 
galaxies in the Universe they allow detailed studies over a 
large cosmic volume. This is highly desirable for a cosmo- 
logical study given that it diminishes the effect of sample 



2.1 Redshift Selection 

The redshift estimates for this above sample were con- 



structed by using the redshift output as given by ANNz ( Col- 



lister & Lahav 2004) an Artificial Neural Network code. This 



empirical photometric redshift estimator learns an effective 
parameterisation of redshift with varying galaxy magnitudes 
(here in u, g, r, i and z bands) by working on a represen- 
tative training set. For our training set we use a subset of 
the ~ 13, 000 spectroscopic redshifts from the 2dF-SDSS 
LRG and Quasar survey (2SLAQ - Cannon et al. 20061, a 
S w ° (declination) stripe within the DR7 imaging area. 
Specifically, we use 5, 482 of these objects as the training set 
with the rest utilised for testing. For this reason and for this 
specific galaxy sample over the range of redshifts of interest 
(0.45 < 2 < 0.65) |Abdalla et al.|2008| found the ANNz train- 
ing method to have the best performance on an evaluation 
LRG sample compared with other redshift estimation codes, 
with average scatter a z = 0.0575 and <r z defined by, 



Uz ~< (Zpho 



>'■ 



(1) 



The performance of this procedure and the representative 
photometric-spectroscopic scatter can be seen in Figure 2 
of |Blake et aLl[2007l and Figures 2-5 in |Abdalla et~aLl[2008l 
The reliability of the neural network training procedure de- 
pends on the training set being completely representative 
of the target galaxy sample. It is noted that by applying 
this 2SLAQ stripe to the wider photometric LRGs there is 
an extrapolation of the training set with sky position. The 



discussion of this potential systematic, however, is left to 
Section 16.31 



2.2 The Colour Selection 

Our SDSS pre-selection and secondary colour selection of 
galaxies is based on and described in Collister et al.|[2007| 
|Blake et al.||2007l and |Cannon et a l. 2006 For example, at 
the start of the 2SLAQ survey there was an alteration in 
the selection criteria used to extract the homogeneous LRG 
sample from the overall galaxy and object population. This 
is associated with the de Vaucouleurs model magnitude id e v 
and also d pcrp , a colour cut which is related to the g, r and 
i model magnitudes via, 



(r-i)-(g-r)/8.0. 



(2) 



We prefer to act cautiously in order to analyse a galaxy 
sample that most represents the training set used to infer 
its properties. Therefore, we use the colour cuts id e v ^ 19.8 
and dporp 0.55 to select and extract the LRG population 
given that these were the selection criteria used in the strict 
majority of 2SLAQ. Again, these cuts were also used in the 
earlier MegaZ-LRG analysis ( |Blake et al.|2007 l. 



2.3 M-star Contamination 

Comparison with the spectroscopic 2SLAQ survey verifies 
that the pre-selection of SDSS galaxies and further colour 
cuts from ICollister et al. 1 120071 and IBlake et all 120071 are 
accurate. Star-galaxy separation was also ensured using the 
criteria ( |Collister et al.||2007[ ), 



'psf ^modcl 



> 0.2 (21.0 - idev) 



-band de Vaucouleurs radius > 0.2 ' 



(3) 



(4) 



However, the presence of M-stars still persist and represent 
the main source of object contamination (ss 5%) within the 
remaining sample owing to similar broadband colours. Gen- 
erally, an uncorrelated sample of stars will act to suppress 



the power of fluctuations (Huterer et al. 20011. However 



one would expect a slightly correlated variation of stellar 
material through the Galactic plane and hence our survey 
area. We therefore remove a large proportion of these con- 
taminants with an extra cut in star-galaxy separation. The 
ANNz code has a star-galaxy parameter 5 sg as an additional 



optional output ( Collister fc Lahav|2004 1. This uses a vari- 
ety of further inputs, including angular size and gauges of 
the light profile throughout the object. These extra parame- 
ters are detailed in Figure 3 of |Collister et al.|2007| The 8 sg 
output parameter varies continuously from 'guaranteed' star 
S sg = to 'certain' galaxy S sg = 1. We remove all objects 
with 5 sg < 0.2, in the processes decreasing the contamina- 
tion fraction to w 1.5% with minimum loss of real LRGs 
( Collister et al.|2007| and |Blake et al.||2007| ). 



2.4 The Angular Selection Function 

The angular selection function, which is used to deter- 
mine the observable boundaries of the survey, was obtained 
from tsChunk . dr7 . best . par downloaded at www . sdss . org/ 
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Figure 1. The SDSS Data Release 7 (DR7) photometric LRG coverage. At 7746 deg 2 it covers 723, 556 galaxies over a redshift 
0.4 < z < 0.7. The three excluded stripes (76, 82 and 86) are visible towards the boundary of the plot. The 2dF SDSS LRG and Quasar 
(2SLAQ) survey and training set constitutes a narrow stripe (i Ri 0°) that passes approximately through the middle of the equatorial 
coordinate system and the bottom of the defined survey. The fluctuations in Galactic extinction are also shown across the survey area. 
The magnitude values are represented by dark blue (0.0 — 0.05 mag), light blue (0.05 — 0.1 mag), green (0.1 — 0.15 mag) and red (> 0.15 
mag). The dust is particularly abundant near the edges of the survey indicating the outer boundaries of the Galaxy. We test the effects 
of this extinction by removing regions with > 0.1 mag (15% of the area) in Section |6.2| 



dr7/coverage. We converted the provided great circle coor- 
dinates and the survey's stripe numbers to declination 
and right ascension before undergoing a HEALPix pixeli- 
sation on a sphere ( Gorski et al.|[2005 l. We used a total of 
3,145,728 pixels (12xn a id e xn a ide where n s id c = 512) over the 
entire sky, placing a zero in pixels corresponding to holes, 
gaps or regions not surveyed and a one in genuinely sur- 
veyed pixels. This discrete survey mask was then overlaid 
with the aforementioned LRG catalogue to leave the final 
galaxy map. We further tested this with n s id c = 1024 to ex- 
amine the effects of a pixelised space. After appropriately 
adjusting the estimated Ct (found in Section [3|, by divid- 
ing by the square of the HEALPix window function tuf , the 
pixelisation effect was found to be negligible. 



We imposed an additional constraint on the mask/map 
by excluding the survey stripes 76, 82 and 86, which are 
widely separated from the rest of the contiguous region. 
These segments act to increase the complexity of the sur- 
vey window function and contribute relatively little extra 
galaxies. The resulting survey used for the primary angular 
power spectrum analysis spans 7746 deg 2 and 723, 556 galax- 
ies over a redshift 0.4 < z < 0.7. This is a 30% larger area 
for analysis than the first and previous MegaZ-LRG survey 
( |Blake et al]|2007| and |Collister et al.||2007[ ). Likewise, it is 
significantly more expansive than the earlier Padman abhan| 



et al. 



2005 



which covering 3,528 deg 2 and 0.2 < z < 0.6 
represents a slightly different LRG population and analysis 
method. The final sky coverage is shown in Figure [I] 



3 THE POWER SPECTRUM MEASUREMENT 

The measurement of the angular power spectrum is per- 



formed by undertaking a spherical harmonic analysis (Pee- 
bles|1973[ ). By explicitly summing the discrete galaxies over 



the incomplete sky we follow the derivation, methodology 
and/or notatio n of |Peebles||1973| |Wright et al]|1994| |Blake 
let al. 120041 and IBlake et al.1l2007l 



One connects the underlying density field in a redshift 
band to the relevant statistical entities by first projecting 
the mass distribution cr(6,(f>). This distribution is then de- 
composed into a series of spherical harmonics Y;, m and their 
corresponding coefficients ai im , 



(5) 



;=o m=—l 



The statistical distribution-the angular power spectrum 
C(—is then given by the multi-realisation expectation of 
these a;, m coefficients, < |a;. m | 2 >. For a full sky survey 
these coefficients represent an orthogonal and normalised 
basis and are thus found by a summation of the spherical 
harmonic conjugate over the galaxy catalogue, 



(6) 



However, in reality one will observe a masked and there- 
fore incomplete sky. This effectively correlates the spherical 
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harmonic coefficients and induces the correction and adjust- 
ment for ioss of power given by, 



\A, 



Jl 



AO, 



(7) 



where N is the number of galaxies, Af2 is the area of the 
sky and the Ii >m and J; >m integrals in Equations [8] and [9] 
are evaluated over the geometry of the discrete survey area. 
I.e., 8Q = 1 for a surveyed pixel and 8fl — for an unsur- 
veyed pixel. The last subtracted term is a correction for the 
statistical distribution of shot noise and is equivalent to the 
expectation of the corresponding harmonic coefficient for a 
random unclustered sample. 



Y,* dfi 



Jl.m — 



IV, 



(8) 



(9) 



One can then obtain the resulting angular power spectrum 
for a given multipole I via an averaging of Ci, m over the 
(21 + 1) di, m values, 



Cobs 
£ — 



psky 



-I °i,m 



21 + 1 



(10) 



The angular power spectrum is independent of m for sta- 
tistical isotropy. The Ce values are further averaged into 
bins of width Al = 10. As seen later this has the effect of 
decorrelating the measurements and provides a more Gaus- 
sian likelihood. We weight this average by the corresponding 
number of ai m s, 



j: e e : +A \2£- 



l)Q obs 



E,; +A ^+i) 



(ii) 



The angular power spectrum in these Al bands is measured 
up to I — 500. One can therefore use these statistics for 
each redshift band within the survey volume. We measure 
the clustering distribution in four such photometric redshift 
bins, each having width Az — 0.05 from z = 0.45 to z — 
0.65. It is this expression in Equation 1 1 1 1 that is presented 
in Fi gure [4] and S ection |3.2| These procedures are in line 
with |Blake et~a l. 2007 and therefore a direct MegaZ-LRG 
consistency check can be made. 

The aforementioned redshift bins are correlated, how- 
ever, as photometric errors scatter galaxies between the bins. 
A small modification to the angular power spectra, 



Cl' J = 



21 + 1 



(12) 



enables a measurement where the harmonic coefficients in 
bin i and bin j have been adjusted for incomplete sky cover- 
age as detailed above. The results are listed in Section [3. 2| 
Note there exist other analogous procedures for 
the analysis of galaxy clustering including, for example, 
quadratic estimators, maximum likelihood methods and ex- 
plicit reconstructions of the power spectrum (E.g. Huterer| 



2003, 



et al.||2001| |Tegmark et al.|[2002| |Seo fc Eisenstein 
Tegmark et al.||2004[ |Blake fe Bridle||2005| |Tegmark"et al 
2006[JPadmanabhan et al.|2007||Blake et al.|2007| and |Re'id 
et al. |2009l. 
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Figure 2. The averaged reconstruction of the input C( field for 
1000 simulations. The solid lines represent the input cosmology 
for the four equally spaced redshift bins (Az = 0.05) between 
z = 0.45 (top profile) and z = 0.65 (bottom profile). The triangles 
are the binned measured values from Equation |11| The plot has 
been truncated at I = 160 as a visual aid; the behaviour beyond 
this point continues in an identical fashion. The accuracy and 
consistency of the code and measurement procedure is clear. 



3.1 Simulations and Gaussian error 

The methodology described above, for the measurement of 
the angular power spectra, was applied to simulated data 
in order to test the procedure and the code. This was per- 
formed by first constructing a Gaussian random field for 
some input cosmology, using best fit WMAP parameters 



(Larson et al. 20101, and subsequently reconstructing this 



cosmology for each of the four galaxy clustering redshift 
bins to be measured. We randomly selected the full set of 
spherical harmonic coefficients a t i m from Gaussian distri- 
butions with widths given by the underlying known cosmol- 
ogy [(CV) 5 ]. The relation between the underlying matter 
power spectrum and the theoretical angular power spectrum 
is described i n Section [4| T hen, using the HEALPix func- 
tion alm2map jGorski et al. 2005) we simulated a pixelised 
galaxy map from these quantities and sampled objects as 
a Poisson realisation of the field. The full angular selection 
function of the survey (Section 2.4 Figure [T]) was imposed 
on the simulated map and the number of galaxies sampled 
in each bin were matched to those present in the observed 
catalogue. This mock data was then analysed with the mea- 
surement pipeline in the same manner as the real data and 
averaged over 1000 simulated realisations. The accuracy and 
reliability of the code and the power spectrum measurement 
procedure is evident in Figure [2] 

One can also use these simulations to check the ana- 
lytic statistical error in the galaxy clustering measurements 
a(Ci), which we assume to be Gaussian. This is extracted 
from the standard deviation over the 1000 mock realisations 
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Figure 3. Left Panel: The analytic Gaussian expression (Equation |13| dashed line) is accurately traced by the 1000 realisation simulated 
error in a redshift band (solid line), shown here for bin 1 (0.45 ^ z ^ 0.5). This demonstrates the approximate validity of the Gaussian 
expression. Right Panel: The agreement is further highlighted by the ratio of the analytic and numerical estimations of the statistical 
error, where the overall behaviour is consistent with unity. The two panels are shown for the first bin only but are representative of all 
other bin combinations. 



at each £. The analytic expression (E.g. Dodelson|2003 and 
Blake et al!]|2007[ ) is given by, 



*(c,; 



/sky (2/ 



(13) 



where f s k y is the fraction of sky surveyed, Afl is the area, N 
is the measured number of galaxies in the bin and Ce is the 
observed or theoretical angular power spectrum. The first 
and second terms in Equation [13] include the necessary er- 
ror contributions from both cosmic variance and shot noise, 
respectively. It also accounts for the reduced error given the 
combination of 2£+ 1 CV, m values into the determination of 
each Ci . For the statistical error in the cross power spectrum 
this generalises to, 



/sky(2£+l) 



( cl + n~Jaq) ( 



ci + 



i 



(14) 

We find the expression reconstructs the simulated error ac- 
curately in each of the four redshift bins and across the entire 
range of £. This is easily seen in the left panel of Figure [3] 
The error ratio, typified by the first redshift bin, is displayed 
in the right panel. Therefore, for our cosmological analyses 
we use the analytic expression only in the covariance matrix, 
evaluated with the model Ces. 



3.2 Results 

We have constructed the galaxy clustering angular power 
spectra Ce for SDSS MegaZ-LRG (DR7), an extension to 
the earlier analysis ( Blake et al.|2007 1 of the original MegaZ- 
LRG catalogue ( |Collister et al.||2007| ). Including 723,556 
photometrically determined LRGs and encapsulating 7746 



deg 2 the measured values in four redshift bins extending 
Az = 0.05 in redshift, from 0.45 to 0.65, are included on- 
lintQ These values are illustrated in Figure [ZlThe full mea- 
surement procedure was detailed in Section[3] In addition, 
the measured cross power spectra between bins are described 
in Section [6T] and are also included online. 

With the angular power spectra the statistical errors 
<j(Ci) on each power spectrum measurement are also in- 
cluded as given by Equation|13| but calculated with the mea- 
sured Ci. They have been further weighted over the Al —10 
band. This was shown in Figure [3] to be a good approxima- 
tion. Note that for the cosmological parameter estimation in 
Section [5] we utilise the Gaussian expression but evaluated 
with model Ces. 

In addition to the simulations described previously we 
also test the measurement pipeline by reconstructing the 
observed Ci as found in the DR4 catalogue. We find these 
values to be identical to lBlake et al. 1 12 0071 As hinted in the 
DR4 results we find that DR7 also exhibits an excess of 
power over the largest scale (l ~ 6 band) in the furthest 
redshift bin. The effects and potential cause of this will be 
discussed in a companion paper. 



4 THEORETICAL POWER SPECTRUM 

In order to deduce the cosmology to match the measured 
angular distribution above one must first have a method for 
connecting the underlying 3D mass distribution to Ce- The 
outline description below simply follows the approach and 



1 http://zuserver2.star.ucl.ac.uk/~sat/MegaZ/MegaZDR7.tar.gz 
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Figure 4. The measured (auto) Angular Power Spectra (Ci) for the photometric SDSS MegaZ-LRG (DR7) population evaluated 
using Equation [Tl] The error bars correspond to those calculated with Equation |13| using the measured power spectrum. These include 
contributions from cosmic variance and shot noise, while accounting for the fraction of the sky surveyed. The solid line is evaluated for 
the the best fit parameters found in Section [5] using the |Smith et al.|2003l non-linear prescription. The panels are: Bin 1 (top left), Bin 
2 (top right), Bin 3 (bottom left) and Bin 4 (bottom right), containing 259,498; 237,564; 155,293 and 71,201 galaxies, respectively. In 
the furthest redshift bin an excess of power is observed over the largest scale. 



notation of Hu terer et al.|[200T| [Tegmark et al.|20 02, Blake 
|et al.|2007| and, most clearly, |Padm^nabharreTal*T2*007 

One starts by noting that before the statistical decom- 
position of the density field into spherical harmonics in Sec- 
tion[3]the field was projected. The same procedure is initially 
followed for the theoretical angular power spectra with the 
3D mass distribution S projected along the line-of-sight 5 2D . 
This gives, 

5 2D =i l J -0^5{\L)W L {k), (15) 

where 8 has also undergone a Fourier transformation. The 
resulting spherical Bessel function j(kz) and the projection's 
weight f(z) have been absorbed into the window function 
given by, 



Wi{k) = J f(z)ji(kz)dz. (16) 

The weight naturally depends on the normalised redshift 
distribution of the objects under consideration f n(z) dz = 1 
and the linear growth factor D(z), 

f(z) = n(z)D(z)(~) (17) 

with the Jacobian relating to the radial comoving coordinate 
x. Using the definition of the power spectrum P(k) for the 
3D density field <5(k), 

< <5(k)<f (k') >= (27r) 3 5 3 (k - k')P(fc) (18) 

the angular power spectrum Ci is found and similarly de- 
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fined to be, 



=< S 2D 6* 2D >-- 



4tt / A 2 (fc)Wf(fc) 



(1A: 



(19) 



The spectrum has been recast into the dimensionless power 
spectrum defined in Equation |20| This power spectrum de- 
scribes the variance of the matter field in logarithmic bands 
and so the equation for Ci is subsequently a weighted in- 
tegral of this quantity over logarithmic intervals (dk/k — 
dlnfc). 



A 2 (fc) 



47rfc 3 P(fc) 
(2tt) 3 



(20) 



This can be further written in terms of the galaxy power 
spectrum with the addition of a linear galaxy bias b, 



P g (k) = b 2 P(k). 



(21) 



For an analysis between redshift bins the above outline can 
be easily extended. The cross correlation of two distinct pro- 
jected mass distributions < 8 2D 6* 2D > leads simply to a 
slight modification in Equation |19| with the window func- 
tion for each bin treated separately, 



C; j =4tt / A 2 (k)W l (k)W j (k) 



dk 



(22) 



For £ > 60 the exact expression (Equation |19[ ) can be sim- 
plified by the small angle approximation (e.g. Blake et al. 
2007b, 



P(k,z) 



n(z) 



/ dx\ - 1 
\~dz) 



dz. 



(23) 



On larger scales (smaller I) this approximation becomes in- 
valid as it seriously underestimates the power in Ce and is 
not used in this regime. In fact, even the exact expression 
does not capture the shape of the true power spectrum be- 
low £ ~ 60. The main reason is because of redshift space 
distortions, which lead to a significant boost in the angular 
power spectrum. 



4.1 Redshift Space Distortions 

The peculiar velocity of a galaxy will cause it to appear 
shifted along the line-of-sight in redshift coordinates (E.g. 

[19801 |Kaiser||1987[ |Fisher| 



Sargent fc Turner||1977l |Peebles | 

et al.||1994| |Heavens fc Taylor||l995f |Hamilton||1998| and 



Guzzo et al. 12008 1. This is relative to the same galaxy car- 



ried along only by the background Hubble flow. That is, the 
redshift distance s of a body will be altered from its true 
distance r, by its own peculiar velocity v = f.v, radially 
from the observer, 



s = r + r.v = r + v. 



(24) 



In redshift space this deviation alters the apparent clustering 
of galaxies and collectively the effect is said to be the result 
of redshift space distortions. 

To include redshift space distortions in the angular 
power spectrum the window function We (k) in Equa tion [T9~ 
is modified suc h that Wi{k) -» We(k) + W^ jk) (E.g. [Fisher 
et al.||1994 and Padmanabhan et aL]|2007 l. This is a result 
of writing the weight properly as a function of redshift dis- 
tance /(s) and assuming that the magnitude of the peculiar 
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Figure 5. A range of theoretical angular power spectra for the 
lowest redshift bin used in this survey (0.45 ^ z ^ 0.5). This in- 
cludes the small angle approximation (Equation |23| dotted line), 
the exact expression with no redshift space distortions (Equa- 
tion [l9] dot-dashed line), the exact expression including redshift 
space distortions (Equation |27| solid line) and also with the ad- 
dition of the partial sky mixing matrix convolution (Section |4.2| 
Equation |28| dashed line). The input parameters are taken to be: 
tt b = 0.05, Q m = 0.3, h = 0.75, <r 8 = 0.8 and b = 1 for all 
four profiles. The small angle approximation is used for multipole 
scales £ > 60 for faster computation in the cosmological analyses. 



velocities are small. This is because with this assumption 
one can perform a Taylor expansion of the weight, 



/( 8 )«/(r) + -^(v(rf).f) 



(25) 



The subsequent window function (remembering Equa- 
tion [l6l therefore now has the two components, We(k) + 



Wp(k), with the latter currently a function of v from above. 
The Fourier transform of v is in turn related to the density 
perturbation through the linear continuity equation, 



v(k) = -i/35 g (kV 



(26) 



with the constant of proportionality (3 known as the red- 
shift distortion parameter. This is commonly approximated 
by p w O.Jn/b, with 7 = 0.55 in LCDM. Substituting this 
into the expression for the window function and Legendre 
transforming (see Padmanabhan et al.||200T for further de- 
tails) eventually leaves one with, 



W t R (k)=/3 J f(y) 
(21 - l)(2l + 1) 



(2Z + 3)(2Z- l) JlK y ' 
3i-2{ky) ■ 



(2« + l)(2Z + 3) 



jl+2(ky) \ dy. 



(27) 



For large values of £ the integral within Equation [27] 
tends to zero and so the total window function is reduced 
to the previous form. In this way, even with the inclusion 
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Figure 6. A slice through the mixing matrix R t t i is plotted 
for two fixed multipole values given by £' = 200 (solid curve) 
and £' = 260 (dashed curve). The amplitude of the matrix peaks 
at those fixed values and decays rapidly within the size of a A£ 
band. This establishes how little correlation is induced by the 
survey's window function. Furthermore, the behaviour is observed 
similarly across all angular scales. Note that the matrix profiles 
have been normalised to unity at their peaks and the vertical axis 
is in logarithmic space. 



of redshift distortions, the small angle approximation is an 
efficient and accurate estimate of the angular power spec- 
trum at small scales. The behaviour of this approximation 
and the effects of the redshift space distortions on the an- 
gular power spectra are illustrated further in Figure [5] In 
addition, we later recast our constraints into limits on the 
distortion parameter f3 as can be seen in Section |5.4| and 
Figure [12] 



4.2 The Mixing Matrix: Partial Sky Convolution 

An additional alteration in the shape of Ce at large scales is 
to account for the partial sky coverage of the real survey. As 
stated in Section [3] this correlates the usually orthonormal 
spherical harmonic coefficients, effectively creating a depen- 
dency on neighbouring scales. The net effect is to slightly 
suppress the shape of the power spectrum Ce below I ~ 60 
as seen m Figure [5] The effect can be calculated by convolv- 
ing with the mixing matrix R^i 



jHauser fc Peebles||1973 
Hivon et al.|2002| and [Blake et al.|2007| ), 



(28) 



The mixing matrix can be pre-calculated and depends purely 
on the survey geometry. It is described by, 
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Figure 7. The spectroscopic redshift distribution n(z) for each 
photometric bin in DR7 is illustrated as a series of histograms. 
Each redshift distribution is fit by a Gaussian function exp[— (z — 
/i) 2 /2<r 2 ], where and a are specified in Tablcjl] The associated 
Gaussian fits are represented by the smooth curves. 



with Wi, the power spectrum of the survey's mask, calcu- 
lated using Equation |30| The 2x3 matrix within R t>l i is 
a Wigner coefficient. For a full sky survey the convolution 
should have no effect on the angular power spectrum and ac- 
cordingly the mixing matrix reduces to the identity matrix 
Rt.t' —> Su' ■ 



Wi 



<r. 



21 + 1 



(30) 



For the DR7 survey geometry the mixing matrix at a given 
£ is seen to be heavily peaked as a function of £' about that 
multipole value. The profile rapidly falls within the chosen 
A£ — 10 bin, implying that only a small correlation between 
the £ bands is introduced by the mask. This is illustrated in 
Figure [6] for two different multipole scales. 



5 THE COSMOLOGICAL ANALYSIS 

We calculate P{k) for the angular power spectrum Ce vls- 
ing camb (|Lewis et al.|2000 i. The halofit fitting function 
( |Smith et al.||2003[ ) is then used to map the linear power 
spectrum into the non-linear regime (large £). To increase 
the speed of calculation we use the small angle approxi- 

> 60 and the full and exact 



mation (Equation 23 1 for 



window function, including redshift distortions (Equation 19 



(29) 



and Equation 27 1 , otherwise. This is all convolved with the 
mixing matrix Re^i as described in the previous subsection. 



5.1 The Redshift Distribution 

The model redshift distribution n(z) in each redshift slice is 
taken to be the form of the spectroscopic 2SLAQ evaluation 
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I" 


a 


Redshift Bin 


Photometric Range 


0.474 


0.0312 


Bin 1 


0.45 < z < 0.50 


0.523 


0.0428 


Bin 2 


0.50 < z < 0.55 


0.568 


0.0433 


Bin 3 


0.55 < z < 0.60 


0.624 


0.0568 


Bin 4 


0.60 < 2 < 0.65 



Table 1. The mean fi and deviation a of the Gaussian fitting to 
the spectroscopic redshift distribution n(z) in each photometric 
bin. This is highlighted in Figure [7] 



set, with the same LRG selection criteria, in that photo- 
metric bin. This is possible because the 2SLAQ evaluation 
objects have both a spectroscopic and photometric redshift. 
These n(z) were fit with a Gaussian function given by, 



n{z) oc exp|^ j 



(31) 



For the cosmological analyses \i and a are fixed to their 
best fit values in each bin. We address this assumption as a 
potential calibration systematic in Section [6. 1| The best fit 
quantities are summarised in Table [I] for the current (DR7) 
data release. In addition, the Gaussian fits to the spectro- 
scopic distributions are illustrated in Figure [7] The vertical 
axis represents the number of spectroscopic 2SLAQ objects 
within a small histogram band (5z). 



5.2 Parameter Constraints: The Single Redshift 
Bins 

We start by undertaking a preliminary cosmological analy- 
sis in each of the four separate redshift bins described pre- 
viously. A conservative choice of parameters is studied such 
that we can test for consistency against the previous MegaZ 

We therefore vary four 



2007) 



LRG analysis (Blake et al. 
quantities: /t, = flj,/fi m , Sl m , as and b; the baryon-to- matter 
density ratio, the matter density, the normalisation of the 
power spectrum and the galaxy bias, respectively. The bias 
is assumed to be scale independent. Along with the earlier 
MegaZ paper the Hubble constant is fixed to Ho = 75 km 
s _1 Mpc -1 and the spectral index to n s = 1. Both as and 
the bias control the amplitude of the power spectrum and 
are thus degenerate with one another. A flat prior is there- 
fore enforced on the former such that 0.7 ^ as ^ 1.1. The 
Universe is assumed to be flat throughout with the equation 
of state fixed to w = — 1. We use all the multipole values 
up to t = 300. This is the scale at which the non-linear cor- 
rections become increasingly significant. For the parameter 
exploration we use the publicly available CosmoMC pack- 
age ( |Lewis fc Bridle|2002| >. 



5.2.1 Data Release 4 

To test for consistency we first perform the cosmological 
analysis on the previous DR4 angular power spectra found 
in |Blake et al.|2007| We find a remarkably similar agreement 
to the previous study in the first three redshift bins over a 
redshift range 0.45 ^ z ^ 0.6. However, for the final and 
furthest redshift bin (0.6 < z ^ 0.65) a large discrepancy is 
discovered when all angular scales to / max = 300 are utilised. 
It is interesting that for this particular redshift bin a large 
excess of power is observed in the measurement of the Ce on 



the largest angular scale (£ ~ 6 band). Even though this is 
approximately at the turnover scale of the power spectrum, 
where one might expect the power to start decreasing, the 
excess was found to be over la from the best fit Ce pro- 
file. One might not therefore expect this anomalous point 
to cause any significant alteration in the cosmological anal- 
ysis. It is important to remember, however, that the error 
on this data point, assigned in the previous study, was the 



error given by the Gaussian expression (Equation 13 1 using 
the data value for the Ce- As the magnitude of this point 
is so much larger than the Ce corresponding to a smooth 
fit through the other data points, the associated data error 
bar is made to appear much larger also. In the parameter 
estimation performed here and in Blake et al. 2007 the er- 
ror and therefore covariance matrix are evaluated using the 
model errors. This is because in a Bayesian analysis one im- 
plicitly assumes the model to be true. Any model spectrum 
attempting to fit the other data points will assign a theo- 
retical value at the largest angular scales much lower than 
that measured and subsequently the error bar will be much 
smaller. Therefore we find the excess power a much poorer 
fit than was ascribed previously in |Blake et al.|[2007| In or- 
der to try and replicate the original DR4 constraint for this 
furthest bin we remove this irregular point. In addition, we 
also follow an analysis using the data errors in the covari- 
ance matrix while including the excess power quantity. The 
resulting contours are shown in the left panel of Figure [8] 

The plot highlights that the excess power at low mul- 
tipoles is indeed significant, with the inclusion of the low- 
est point dragging the constraint to much lower values of 
Q m (red contour). Also, the figure reiterates the notion that 
the data error (green contour) acts to buffer against this 
anomaly given that the contour is similar to the model anal- 
ysis that excludes the excess power (blue contour). When 
fitting with the model errors and no excess power we find 
the constraints to be identical to those in [Blake ct al 1 120071 
and also consistent with the three other redshift bins. 



5.2.2 Data Release 7 

The angular power spectrum for the last redshift bin was 
measured for DR7 in Section |3.2| and shown in Figure [4] 
Once again an excess of power is detected at this high red- 
shift. However, there seems to be a slight hint of an ease in 
tension as the magnitude of the DR4 point is found to be 
40% higher than the newly measured DR7 value. We there- 
fore undertake a cosmological run for this bin using the ex- 
cess power point and also with it removed to test the effects. 
We find that despite the more recent decrement in the Ce 
on these large scales the inclusion of the quantity still sig- 
nificantly affects the parameter constraints found with the 
bin. This is illustrated clearly in the right panel of Figure [8] 
Again, with this point excluded the fourth redshift bin is 
found to be consistent with the other three slices. For the 
new DR7 release the associated constraints for every redshift 
bin are displayed in Figure [9] and Figure [To] 

We therefore choose to continue the galaxy clustering 
study by excluding the anomalous excess power in the I ~ 6 
band for the furthest redshift bin. It is intriguing that slight 
hints of excess power have also been seen in Padmanabhan 
|et al.|[2007| and in the study of the maxBCG cluster power 
spectrum by [Huctsi 2009. We discuss this signal and poten- 
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Figure 8. Left Panel: Constraints on the MegaZ LRG (DR4) highest redshift bin (0.6 < z ^ 0.65) using model Cg$ in the error expression 
(red/leftmost contour), data Cg& in the error (green/rightmost contour) and model errors with the lowest multipole removed (blue/central 
contour) . The last analysis gives constraints consistent with the previous |Blake et al. | | |2007[ l study. Right Panel: DR7 constraints on the 
same bin using model errors (red/left contour) and model errors with the lowest multipole band removed (blue/right contour). Despite 
a slight decrease in the excess power in DR7 the observed shift in constraints above show the contribution from the anomalous low band 
to still be significant. The blue contour analysis is consistent with the other three redshift bins (Figures |9| and [Tp| and as such this point 
is removed from all subsequent analyses. 



S7 1) 



Redshift Slice 



0.166 ±0.066 
0.136 ±0.069 
0.206 ± 0.062 
0.146 ±0.076 
0.173 ± 0.0462 



0.253 ±0.049 
0.251 ± 0.051 
0.274 ± 0.052 
0.248 ±0.067 
0.260 ±0.0351 



Bin 1 
Bin 2 
Bin 3 
Bin 4 
All bins 



0.163 ±0.0480 0.234 ±0.0309 



All bins 



Table 2. The marginalised mean values obtained from the anal- 
yses of the galaxy clustering angular power spectra Ci. fb = 
£lb/£lm, f^m, eg and b are varied for each single bin run. In the 'all 
bins' analyses all the bins were combined together using the full 
covariancc matrix and a bias parameter for each bin (&i, 62 > 63, 
64). In the last analysis the lowest multipole band in the highest 
redshift slice is included. 



tially related systematics or causes further in a companion 
paper. 

When obtaining /(,, Q. m or b all the other parameters are 
marginalised over. The bias is subsequently seen to enlarge 
with an increase in redshift. This is partially due to the ob- 
served galaxies in the furthest redshift bin necessarily being 
more luminous, resulting from the pseudo-magnitude limit 
in the survey. They are therefore observed to be more highly 
clustered ( |Blake et al.j2007| . All the inferred constraints are 
summarised in Table [2] 



5.3 Parameter Constraints: The Combined 
Redshift Bins 

We now combine the data from each of the four redshift 



bins. These bins are not independent, however, as photo- 
metric redshift errors act to disperse galaxies throughout 
the bins. Another way of noting this is to observe that the 
Gaussian redshift distributions, as seen in Figure [7] over- 
lap for each bin. We therefore use the full covariance matrix 
in the analysis. The variance element corresponding to the 
same redshift bin (e.g. between C\ and C\) is given by the 
square of Equation [13] using the theoretical expression for Ct 
as before. The covariance elements between different bins are 
described by, 



Cav(Cj,Cf 



/ sky (2£± 



(32) 



In this way the whole matrix allows for the covariance be- 
tween all bin combinations but not multipole bands. This 
is a good approximation given our earlier discussion of the 
highly peaked mixing matrix Rijt (Figure [6j|. 

We include a redshift dependence in the galaxy bias, 
to the extent that each redshift bin is assigned a separate 
bias parameter (61, &2, &3 and 64) in the cosmological run. 
Potentially there could also be added complexity in the bias 
(E.g. |Swanson et~al][2008| and |Cresswell fc Percival|[2009| 



but this is beyond the scope of this current work. 

The marginalised best fit parameters are again listed 
in Table [21 with the corresponding contours displayed in 
Figure |ll[ As found with the individual bins the contours 
can be seen to visibly rise along the bias axis with an increase 
in redshift. Moreover, the four bias quantities are seen to be 
high implying that the LRGs strongly trace the underlying 
mass distribution. In particular we find: fei = 1.47 ± 0.15, 
b 2 = 1.71 ± 0.17, 63 = 1-80 ± 0.18 and 64 = 2.05 ± 0.21. 
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Figure 9. MegaZ LRG DR7 constraints on fi, = Of,/t2 m and f2 m for four separate redshift bins, b and erg have been marginalised over 
and Ho and n s are fixed to 75 km s~ 1 Mpc -1 and 1, respectively. The panels are: Bin 1 (top left), Bin 2 (top right), Bin 3 (bottom left) 
and Bin 4 (bottom right). The inner and outer contours are the 68% and 95% confidence levels, respectively. 



Finally, in the top panel we include a calculation of the 
combined bins with (red/left contour) and without (all blue 
contours) the lowest multipole band measured in redshift 
bin 4. As with the individual bin the excess power is seen 
to systematically displace the marginalised distribution and 
once again is removed from all other constraints. 



5.4 Measuring Redshift Distortions 

We noted previously that the peculiar velocity of a galaxy 
will cause it to appear shifted along the line-of-sight. The 
modification to the window function was described and de- 
tailed in Section [4. 1| For angular power spectra Ce this can 
cause a significant effect over large scales as seen in Fig- 
ure [5] Using this sensitivity we recast the bias parameter(s) 
in our MCMC chains into the redshift distortion parameter 
£(z) ~ Q m (z) 55 /b, where we take z to be the mid-point 
of the bin in question and also b is the bias of that bin. 



For the combined bin analysis this renders four distortion 
parameters ft = £(0.475), £ 2 = £(0.525), /3 3 = £(0.575) 
and £4 = £(0.625). The resulting limits from each separate 
bin and from the combined analysis are listed in Table [3] 
and are illustrated in Figure |12| It is interesting to note 
that one could use a similar methodology to test for devi- 
ations to gravity or dark energy clustering by allowing the 
exponent in the distortion parameter to vary as a free pa- 
rameter, e.g. £ « n m (z) J /b ( |Linder|200"7| |Guzzo et al.|2008 



and Thomas et al. 



120091). Further, with this data set and 
method one could also fit £ — (l/b)dln_D/ dlna and/or 7 
independently. 



5.5 Other Studies 

In as much as other analyses can be compared, with varying 
parameter choices and assumptions, these results are concor- 
dant but competitive with respect to recent studies of SDSS 
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Figure 10. MegaZ LRG DR7 constraints on f2 m and the bias b for four separate redshift bins. = Q^/Qm and erg have been marginalised 
over and Hq and n a are fixed to 75 km s^ 1 Mpc -1 and 1, respectively. The panels are: Bin 1 (top left). Bin 2 (top right), Bin 3 (bottom 
left) and Bin 4 (bottom right). The inner and outer contours are the 68% and 95% confidence levels, respectively. 



Separate Bins 


Combined Bins 


Redshift 


0.423 ±0.058 
0.392 ±0.057 
0.440 ± 0.045 
0.345 ± 0.055 


0.475 ± 0.050 
0.418 ±0.043 
0.409 ± 0.042 
0.367 ± 0.038 


0(0.475) 
0(0.525) 
0(0.575) 
0(0.625) 



Table 3. The marginalised mean redshift distortion parame- 
ters measured from the analyses of the galaxy clustering angular 
power spectra Cg . /(, = Qb/^m> ^m> ""8 have been marginalised 



over with a prior of 0.7 ^ erg ^ 1.1 on erg. 
75 km s" 1 Mpc" 1 and 1, respectively. 



Ho and n a arc fixed to 



galaxy clustering. These often include alternate or earlier 
data sets and at different redshifts. This includes IPadmarH 
abhan et al. 2007} an analogous photometric study to [Blake] 
eFal*pob7[ that instead reconstructs the 3D real space power 
spectrum. Apart from these two works, studies have tended 



to focus on the spectroscopic samples, such as Tegmark et al. 
20041 |Tegmark et al.|[20Tj6] |Cabre fc Gaztafiaga||2009| and 



Sanchez et al 



2009 



and most recently [Reid et al.|20 09 (DR7). 
Furthermore, the SDSS galaxies have permitted mea- 
surements of the Baryon Acoustic Oscillations with Percival| 
|et al.|2007||Gaztanaga et al.|2008| and |Percival et al.|2009| 



6 SYSTEMATICS AND FURTHER TESTS 



The earlier MegaZ release ( [Blake et al.| [20?J7[ > performed 
a series of systematic tests based naturally on examining 
variations across the sky. This included astronomical seeing, 
overlapping survey stripes, regions of low Galactic latitude, 
varying completeness and variations in star-galaxy separa- 
tion. All potential effects were found to have little or no in- 
fluence on the estimated power spectra. The aforementioned 
paper also highlighted the impact of photometric errors for 
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Figure 11. MegaZ LRG DR7 cosmological constraints given by the combination of four redshift bins between 0.45 < z < 0.65. The 
red/left 2D distribution in the top panel shows the systematic shift induced by including the excess power measured over large scales 
in the highest redshift bin, whereas the normal blue contours have this anomalous point removed. The bottom four panels show the 
determination of the bias where b\, b2, 63 and &4 are the quantities in sequentially higher redshift bins. These correspond to marginalised 
limits of 61 = 1.47 ± 0.15, b 2 = 1.71 ± 0.17, b 3 = 1.80 ± 0.18 and 64 = 2.05 ± 0.21. 



14 Thomas et al. 



0.55 




°fo5 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 °fo5 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.4 5 0.50 0.5 5 0.60 

Q m Q m Redshift (z) 

Figure 12. Limits on the redshift distortion parameter f) f2 m (z)°- 55 /fe and f2 m for an analysis of each separate redshift bin (left 
panel); for f3\, 02> /?3, Pi and f2 m for the four redshift slices combined (middle panel); and the corresponding marginalised values and 
error derived from this combined bins analysis (right panel). Again, /), and erg have been marginalised over; there is a 0.7 ^ erg ^ 1.1 
prior on o-g, and Hq and n s are fixed to 75 km s — 1 Mpc~ 1 and 1, respectively. 



LRGs given their location on the galaxy luminosity function. 
This function <p(M) describes the number of galaxies that 
have absolute magnitudes M within an interval M + AM . 
The position of the galaxy sample under consideration is 
one where the gradient of this function is high. Therefore, 
any slight systematic shift in M will impart a large system- 
atic shift in the number of galaxies. If this systematic shift 
were some function of sky position, for example, it could 
contribute significantly (and artificially) to the galaxy clus- 
tering signal at that scale. With these considerations in mind 
we therefore choose to examine the redshift distribution, the 
role of extinction and the process of photometric redshift es- 
timation. 



6.1 Redshift Bin Cross Correlations 

A useful test of any known or unknown systematic present in 
the study is the cross angular power spectra (Equation |12| 
and Equation 22 1. A signal in these quantities should be 



the result of photometric errors scattering galaxies between 
bins as predicted by the spectroscopic redshift distribution 
defined earlier and the best fit auto-power spectra. Any sig- 
nificant alteration in the measurement relative to the the- 
oretical CI' 3 could indicate an additional systematic in the 
photometry, extinction correction or an ill-calibrated red- 
shift distribution, for example. 

We measure the cross power spectra in each of the six 
cross-bin combinations (note that C\' 3 = Cf) in multipole 
bands of A£ = 10 up until Z max = 500 as performed previ- 
ously for the auto power spectra. The observed values are 
included online and are plotted in Figure [13] along with their 
associated error bars. The solid lines in these plots show the 
predicted theoretical spectra using the best fit values from 
the combined bins analysis and the corresponding Gaus- 
sian redshift distributions. For nearby bins there is excel- 
lent agreement in the values. However, the anticipated cross 
spectrum between bin 1 and bin 4 (middle left panel) suffers 
from a lack of amplitude and consequently does not fit the 
mean profile of the data well. This is most likely the result 
of the Gaussian redshift distributions being weak fits to the 



spectroscopic profiles far from the mean of the distribution. 
As can be seen in Figure [7] the Gaussian underestimates the 
number of galaxies far from the bin centre. This will lead 
to an under prediction in the cross term. Less dramatic is 
the apparent marginal overestimation of the cross spectrum 
between bins 2 and 3 (middle right panel). This might be 
the result of the Gaussian smoothing adding slightly more 
galaxies in the overlap region between the redshift slices. 

To test this hypothesis we interpolated the spectro- 
scopic distribution with a spline through the n(z) histogram. 
Then using this more 'realistic' profile we re-evaluated the 
theoretical cross power spectra in the bin. These are shown 
as the dashed lines in the cross spectrum panels. For the 
most physically separated bins (1 and 4) this is seen to give, 
as predicted, a boost in amplitude and a better fit to the 
data points. This could hint that the use of a fixed /x and 
a in the Gaussian is not completely optimal. Potentially, 
one could use a splined distribution in the cosmological pa- 
rameter estimation, however it is important to note that 
this could introduce errors of its own. For example, it might 
propagate inherent fluctuations in the profile that are par- 
ticular to that bin and patch of the sky into the analysis. In 
any case the main effect seen from these figures will be to 
slightly alter the furthest (and least contributing) corners of 
the covariance matrix, whereas the majority are good fits. 
Finally, with the redshift function now more fully tested in 
the cross correlation measurement it is interesting to see 
that in several of the bins correlated with bin 4 there exists 
a slight excess of power. This could point towards a residual 
systematic in the catalogue. 

6.2 Extinction 

Light from galaxies is potentially absorbed, scattered or re- 
emitted by the dust and gas within our own Galaxy. This 
Galactic extinction has the capacity to be one of the dom- 
inant systematics in a galaxy survey such as this. For ex- 
ample, extinction can preferentially absorb light at the blue 
end of a galaxy's spectral energy distribution thus making it 
appear redder and more LRG-like. Alternatively, it can have 
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Figure 13. The measured cross Angular Power Spectra (C» J ) for the photometric SDSS MegaZ-LRG (DR7) population evaluated using 
Equation |12| The error bars correspond to those calculated with Equation |14| using the measured power spectrum. The solid lines are 
evaluated for the the best fit parameters found in Section |5.3| using the the Gaussian redshift distributions. The dashed lines are the 
theoretical power spectra using a spline interpolation of the spectroscopic distribution. The panels are: Bin 1,2 (top left), Bin 1,3 (top 
right), Bin 1,4 (middle left), Bin 2,3 (middle right), Bin 2,4 (bottom left) and Bin 3,4 (bottom right). 
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Figure 14. The exaggerated effect caused by neglecting the ex- 
tinction correction for the i^ cv magnitude cut in bin 1 (dashed 
line). Although this is not used in the study it highlights the range 
of scales that could be affected by any such systematic error in 
the correction. The extinction corrected spectrum is shown by the 
solid points with associated error bars. The solid line is a best fit 
profile for comparison. 



the effect of scattering faint galaxies from the sample. As the 
contribution from our own Galaxy changes as a function of 
position (Figure [l]) this is a cause for concern given that we 
are interested in inferring cosmological quantities through 
statistical variations across the sky. Worse still, it could act 
to further systematically bias our redshift estimates given 
that the ANNz derived galaxy catalogue is a spatial extrap- 
olation of the 2SLAQ training set, which confined to a stripe 
at 8 ~ ° , covers a limited region of Galactic extinction. 

Fortunately detailed maps of Galactic extinction are 
available (Schlegel et al. 1998 Figure [T]) and subsequently 
the u, g, r, i and z bands used are dereddened model mag- 
nitudes, i.e. they are extinction corrected. Figure |14] shows 
the exaggerated effect that is the result of not adjusting 
properly for the presence of dust. In this plot the angular 
power spectrum is evaluated for the catalogue when the idev 
magnitude cut is not extinction corrected. This causes extra 
galaxies to be scattered from the sample in different regions 
of the survey area and a large boost of power is observed. 
Although the values used for our galaxy clustering statistics 
are corrected for extinction it could be that there are errors 
in the correction map. If these errors were related to the 
magnitude of extinction or again varied with position, then 
they too would propagate into the LRG sample. 

To test for extinction correction errors we repeat the 
whole measurement of the angular power spectra with re- 
gions of high extinction removed (> 0.1 mag). This consti- 
tutes a 15% removal of the survey area. The changes in the 
power spectra are plotted in Figure [15] with error bars as 
derived before but for the 'cut' spectra. It is clear that the 
general profiles are not significantly affected. This result is 



consistent with the preliminary examination in Blak e et al.| 
[20071 and lAbdalla et aTll2008l 



6.3 Photometric Redshift Codes 

We chose to derive the redshift estimates for our catalogue 
using ANNz given that it was shown in |Abdalla et al.|2008| 
to have the best performance for this survey. This is par- 
tially due to the presence of the specific 2SLAQ training 
set. We noted that this training set was limited in sky po- 
sition, and therefore extinction, and the application to the 
wider survey was an extrapolation. We test this limitation 
now by comparing the results from different redshift codes. 
Note that this also represents a cosmological comparison of 
the codes with current data. 

lAbdalla et al. 120081 evaluated the SDSS LRG DR6 cat- 
alogue with six photometric codes: ANNz, HyperZ (Bol- 



BRA (Fcldmann et al 



zonella et al. 20001, SDSS ( Padmanabhan et al. 2005), Le 



PHA RE (jllbert et al.||2006|> , BPZ ( |Benitez||2000| and ZE 



2006). These other codes do not re- 



quire representative training set data but instead use a va- 
riety of methods including various LRG templates to ob- 
tain the redshifl]^] As the template based procedures do not 
utilise a spatially confined training set one can argue that 
they are effectively blind to this potential redshift calibra- 
tion systematic. HyperZ was also analysed using two dif- 
ferent sets of templates. This includes observed templates 



provided by Coleman et al. 1980 (CWW) and those syn- 
thetically produced by Bruzual & Chariot 2003 (BC). 

We can use the catalogues presented in |Abdalla et al7| 
|2008| as a test of our work because we compare between codes 
and the DR6 survey area comprises the same contiguous 
region but with only a 1% smaller survey area than in DR7 
(7670.9 deg 2 ). To make a fair like-for-like comparison we 
perform a star galaxy separation based on the ANNz output 
as the other codes do not have this option. 

We plot the overall cosmological constraints for the four 
combined redshift bins from these different catalogues in 
Figure |l6| It is reassuring that there is excellent agreement 
between ANNz and several of the template based codes, 
including SDSS and Le PHARE. Moreover, it is interest- 
ing that where there is some difference it is between tem- 
plate based procedures (ZEBRA, HyperZ BC and HyperZ 
CWW) and not with the independent training set method. 
This clearly indicates consistency for our chosen procedure 
and furthermore shows no degrading effects from a limited 
training set. This is consistent with |Abdalla et al.|2008| and 
all of the contours are consistent within la. However, for 
more statistically discriminating surveys in the future the 
template differences may become important. 

Note that we have accounted for the different redshift 
distributions (fi, a) that each code predicts given the varying 
spectroscopic-photometric relation. These distributions are 
given and are illustrated in Table [AT] and Figure [AT] in the 
Appendix, respectively. 



2 Please see Abdalla et al. 2008 for an overview of these codes and 
templates or the specific references themselves for more detail. 
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Figure 15. The difference between the Angular Power Spectra with regions of high Galactic extinction removed (> 0.1 mag) and with 
the whole survey geometry included as before (i.e., cP artlal _ C7 a11 ). This is to test for possible extinction contamination in the analysis. 
The panels are: Bin 1 (top left), Bin 2 (top right), Bin 3 (bottom left) and Bin 4 (bottom right). There is no observable discrepancy 
between the two calculations. 



7 CONCLUSION 

We have measured and constructed the galaxy catalogue and 
angular power spectra for 723, 556 Luminous Red Galax- 
ies in the Sloan Digital Sky Survey Data Release 7 - called 
MegaZ DR7. This photometric extension to the previous re- 
lease ( Blake et al.|2007 l represents the largest galaxy survey 
to date. Covering 7746 deg 2 over the redshift range 0.45 < 
z < 0.65 we find constraints of /{, = fii,/O m = 0.173 ± 0.046 
and flm = 0.260 ± 0.035. We also use our constraints to 
place limits on the redshift distortion parameter /3 given 
that the Ces are sensitive to distortions over low multipoles. 
We find /3 (z=0 .475) = 0.475±0.050, /3 (z=0 .525) = 0.418±0.043, 
/3 (z =o.575) = 0.409 ± 0.042 and /3 {z=ox>25) = 0.367 ± 0.038. 

It is reassuring that the galaxy clustering results on fb 
and Q m are consistent with the most recent WMAP analysis 
| |Komatsu et al.| 2010). This is a crucial consistency check as 
the two surveys probe vastly contrasting cosmic epochs and 



are subject to different systematics. In addition, it seems 
the photometric approach to modern cosmological surveys 
is justifiable with competitive and concordant results com- 
pared to the spectroscopic SDSS survey ( Reid et al.| [2009). 
This high redshift photometric survey could be combined si- 
multaneously with spectroscopic BAO measurements, such 
as |Percival et al.||2007| with no complex cross-covariance 
(i.e. they are independent). This can be highly complemen- 
tary for the parameter space, particularly in neutrino mass 
measurements as shown in |Thomas et al.||2010| The tight 
constraints on the matter densities show there is now over- 
whelming and precision evidence for some dark energy-like 
component to the cosmos when including photometric data 
from the late-time Universe. A more complete and extended 
set of cosmological implications for this data set and in com- 
bination with other data sets are to be included in a com- 
panion paper, along with the likelihood. 
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m m 

Figure 16. Cosmological constraints on the four combined bins using various photometric codes. This represents a comparison between 
the codes and templates with current data. Left Panel: The training set method with ANNz (thick dashed line) is extremely consistent 
with Le Phare (blue/darker) and SDSS (green/lighter) template procedures and so vindicates the spatial limitation of the 2SLAQ training 
set, including with extinction, for example. Right Panel: The only difference is between a subset of template procedures highlighted by 
HyperZ BC (yellow/lightest), Zebra (red/darker) and HyperZ CWW (grey/darkest). However, all contours still overlap within lcr. Again, 
ANNz is illustrated by the dashed contour for reference. 



We examined the possibility of residual systematics 
in the catalogue from photometric calibration, extinction, 
sky extrapolation and photometric estimation procedure. 
We found no significant alteration in the power spectra or 
combined constraints from these tests. In addition, the lat- 
ter analysis included a comparison of different photometric 
codes and templates. We found that ANNz is highly con- 
sistent with other template based redshift estimation proce- 
dures, such as SDSS and Le PHARE. Small differences were 
seen between some of the template methods but all codes 
tested produced results that were overlapping within la. 
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Figure Al. The spectroscopic redshift distributions n(z) in four photometric redshift bins for different photometric codes and templates. 
The y-axis denotes the number of galaxies in some small redshift interval Sz. This is based on the DR6 catalogues provided by |Abdalla| 
et al. 2008 In order to derive the combined constraint for each code each differing redshift distribution was taken into account. For 
the HyperZ code CWW and BC refer to observed and synthetic templates given by |Coleman et al.|1980| and |Bruzual He Ch ariot 2003 
respectively. For BPZ we use the 'Bayesian' output for the redshift and therefore (99) and (95) refer to cuts of odds > 0.99 and 
odds > 0.95, respectively (see |Abdalla et al.|2008] and references therein for details). 



Code 


Mi 


01 




02 




0"3 




(74 


ANNz 


0.475 


0.0323 


0.523 


0.0428 


0.571 


0.0429 


0.627 


0.0537 


SDSS 


0.499 


0.0431 


0.544 


0.0420 


0.573 


0.0383 


0.618 


0.0421 


Le PHARE 


0.509 


0.0463 


0.557 


0.047 


0.599 


0.0489 


0.632 


0.0572 


ZEBRA 


0.499 


0.0403 


0.549 


0.047 


0.578 


0.069 


0.571 


0.085 
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Table Al. The mean fjbi and deviation Oi of the Gaussian fitting to the spectroscopic redshift distribution n(z) in photometric bin i for 
different codes in DR6. The spectroscopic n(z) are illustrated in Figure [Al] 



